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I . INTIiDDUCTION 


For the past several years semiconductor manufacturers have been concerned 
with iiiperfections that are induced during device fabrication. Considerable 
effort has been directed toward understanding the causes of these dislocations, 
and various processing steps have been identified as major contributors to the 
generation of dislocations. During diffusion thermal gradients sufficient to 
cause the generation of dislocations are produced. Impurity depositions lead to 
concentration gradients which also produce stresses that result in dislocation 
creation. It is further known that the thermal mls-match of the oxide at the oxide- 
silicon interface can lead to the generation of undesirable imperfections in the 
wafer. While these steps have been identified as crucial ones in controlling 
process-induced dislocations, the critical conditions present when the dislocations 
are generated during each step have not yet been determined. To understand the exact 
conditions present when dislocations are generated in any one of these steps 
requires that one be able to separate the various Influencing factors. It is 
necessary then to be able to selectively create the dislocations while maintaining 
control of all the critical parameters. One critical parameter which has been 
largely ignored is the gravitational stress in the wafer which is ever present in 
all earth bound processing experiments. 

The primary objective of this research has been to establish the role played 
by gravitational stresses in process-induced defect generation. 



II. One-g Processing Experiments 


A. Experimental Procedures 

The primary objective of these initial experiments has been to demonstrate 
the influence of gravity on creating dislocations in silicon wafers during thermal 
cycling. Thus^ thermal cycling of silicon wafers were performed in a controlled 
ambient where no impurities were present and oxidation was minimum. Both n and 
p type silicon wafers having a diameter of 1.25 to 1.5 inches with {ill} orientation 
and 8--16 Q-cm resistivity were used. These wafers were positioned either individually 
or in pairs on a specially designed quartz carrier which supported the wafers at 
the bottom edge. The carrier has slots which allows the wafers to be supported 
either vertically or at a 45° cantilevered position. These positions are shown 
in Figure 1 which also indicates the ”top" and ^^bottom” surfaces for the inclined 
wafers. 

The surface dislocation densities were measured quantitatively by the well- 
known Sirtl etch technique. The wafers were immersed in the etch solution from 
two to five minutes. This etching procedure, produced an etch pit at the site of 
each dislocation intersection with the wafer surface. The etch pits were then 
counted under microscope inspection at each of the standard ASTM positions. The 
location of the nine standard ASTM positions is shown in Figure 2 for 1.25” 
wafers. This technique requires that the wafer surface under examination be hi^ly 
polished. Several wafers from each batch of starting material that was used were 
Sirtl etched and all were found to have a negligible backgroxmd density of dislocations. 
All the Wafers used in this study were thermally cycled under varying conditions 
described in the following sections and the process-induced dislocations were revealed 
using the Sirtl etch technique. 




45 CANTILIVERED 


POSITI ON 

I ENTATION 





Location of Count Position on Wafer 

Field 1, 6 : 2,8 mm from edge 

Field 2, 7 : 7.3 mm from edge 

Field 3 : 16 mm from edge 

Field 4, 8 : 24.7 mm from edge 

Field 5, 9 J 29,2 mm from edge 

Field of View: 0.11 mm^ 


Figure 2 
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B. Results on Vertically Standing Wafers 

In order to separate as much as possible the influencing factors referred 
to in the introduction a large number of wafers were thermally cycled in the 
vertical position varying; 1) the insertion — withdrawal rate, 2) the furnace 
flat zone temperature, and 3) the duration of the thermal cycle. Some of these 
results are shown in Table I for 1250‘’C and Table II for 1050°C, These results 
indicate that an insertion - withdrawal rate of up to 2 in/min is sufficiently 
slow to minimize defect generation due to thermal gradient stresses with the 
particular low thermal mass carriers that were used in all these experiments. 

The determination of this maximum insertion - withdrawal rate was one of the 
primary objectives of this series of thermal cycles. However the most significant 
features of the results shown in these Tables is the Increasing damage done as 
the cycle duration is Increased. Also. the amotmt of damage done for a given cycle 
duration increases rapidly with increasing temperature. For example, about the same 
amount of damage is done in 12 to 15 minutes at 1250®C as in one to three hours at 1050"C 

C. Results on Wafers Cycled in the 45* Cantilevered Position 

Using the 2 in/min insertion-withdrawal rate established during the vertical 
cycling experiments, wafers were cycled in the 45° cantilevered position in the 
following ways: 1) Two wafers, one with the polished side up and one with the polished 

side down 2) Wafers with both sides polished. It is Important to note in this 
case that the gravitational stress produced a compresslonal deformation on the wafer’s 
bottom side and a stretching deformation on its top side. These experiments were 
run with furnace temperatures ranging between 1000°C and 1250°C. After thermal 
cycling, each wafer was given a 2 minute etch to remove the (very thin) oxide. Then, 
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Table I 





W A F E R S 

CYCLE 

D V E 

R T I C A 

L L Y AT 

12 5 

p® C 




Sample 

Temp. 

C 

Stay 

(Time) 

Push in 
(Withdraw) 
rate 

Loc. //I 
Etch 
Pits 

Loc. #2 
Etch 
Pits 

Loc. #3 
Etch 
Pits 

Loc.#A 

Etch 

Pits 

Loc . #5 
Etch 
Pits 

Loc. //6 
Etch 
Pits 

Loc. #7 
Etch 
Pits 

Loc. #8 
Etch 
Pits 

Loc.#^ 

Etch 

Pits 

' #39 

1250 

15 sec 

very fast 

10^ 

2 

0 

4 

96 

2 

0 

0 

103 

#40 

1250 

0 

Ain /min 

3 

10 

2 

0 

75 

0 

2 

6 

78 

#41 

1250 

0 

2 in /min 

10 

9 

0. 

0 

83 

0 

2 

12 

97 

#42 

1250 

0 

lin/min 

14 

8 

0 

9 

25 

7 

2 

0 

107 

#43 

1250 

15 sec 

. lin/min 

27 

12 

37 

21 

39 

29 

39 

41 

143 

#44 

1250 

15 sec 

24in/min 

0 

0 

0 

0 

0 

1 

0 

1 

22 

#45 

1250 

1 min 

2 Ain /min 

0 

2 

0 

0 

13 

1 

0 

42 

54 

#46 

1250 

3 min 

2 Ain /min 

10 

24 

0 

15 

10 

2 

0 

0 

37 

#47 

1250 

6 min 

2 Ain /min 

102 

2 

3 

3 

116 

5 

0 

86 

142 

#48 

1250 

12 min 

2 Ain /min 

6 

7 

20 

68 

6 

63 

7 

52 

78 

#49 

1250 

12 min 

2 Ain /min 

23 

6 

1 

7 

33 

15 

5 

41 

4x1 02 

#50 

1250 

15 min 

2 Ain /min 

107 

6x10^ 

4x10^ 

6x10^ 

2x10^ 

4x10^ 

98 

3x10^ 

5x102 

#51 

1250 

15 min 

2 Ain /min 

10 3 

10 3 

5x102 

6x102 

10 3 

63 

92 

10 3 

10- 

#52 

1250 

20 min 

2 Ain /min 

2x1Q2 

5x10^ 

2x1Q2 

2x10^ 

4x10^ 

3x10 3 

3x103 

10 2 

10- 

#54 

1250 

12 min 

2 Ain /min 

32 

5 

60 

73 

92 

7 

21 

5x102 

10- 

#55 

1250 

13 min 

2 Ain /min 

21 

2x10^ 

44 • 

3x10 3 

10 3 

5x10^ 

6x102 

103 

5x10- 

#56 

1250 

14 min 

2 Ain /min 

23 

62 

52 

4x10 3 

3x10 3 

7 

15 

103 

10- 

#57 

Notes; 

1250 

15 min 

2 Ain /min 

13 

3x10^ 

112 

4x10 3 

53 

22 

6 

103 

10- 


(1) The numbers on this table are the etch pits counted within the field of view (O.llinm^), 

(2) All wafers were Sirtl etched for 4 minutes. 

(3) All wafers were vertically placed on low mass quartz boat while baking. 



Table II 


Sample 

Number 

Temp. 

w_ 

A F E R 

S CYC 

LED V 

E R T I C 

ALLY 

AT 10 

5 0° C 

Log. # 7 
Pits 

Log. # 8 
Pits 

Log. # 9 
Pits 

Time 

Stay 

Rush 

Rate 

Loc.#l 

Pits 

Log. # 2 
Pits 

Log. # 3 
Pits 

Log. # 4 
Pits 

Log. # 5 
Pits 

Log. # 6 
Pits 

#58 

1050 

30min 

Imin 

3 

12 

6 

11 

39 

7 

7 

20 

23 

#59 

1050 

Ihr 

Imin 

34 

45 

38 

29 

46 

88 

69 

1Q2 

3x102 

#60 

1050 

Ihr 30min 

Imin 

20 

6 

1 

7 

39 

25 

7 

49 

4x10^ 

#61 

1050 

2hr 

Imin 

1Q2 

60 

5x10^ 

79 

4x10^ 

102 

10^ 

46 

10 3 

#62 

1050 

Ihr 30mln 

Imin 

35 

77 

10^ 

67 

94 

45 

7 

10^ 

10 3 

#63 

1050 

2hr 

Imin 

10^ 

10^ 

57 

9 

6x10^ 

102 

31 

10^ 

1Q2 

#64 

1050 

3hr 

Imin 

10 3 

10 3 

5x10^ 

60 

10^ 

lO^ 

66 

10 3 

3x10^ 

#65 

1050 

4hr 

Imin 

102 

10 3 

10^ 

10^ 

10^ 

10 3 

57 

10 3 

10 3 


View at each point is a circle, 0.375inm in diameter, or Silicon wafer, 
boron doped { 111] orientation, flat at (110)- All wafers Sirtl etched before 
baking. 
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a sequence of Sirtl etches were performed with a total Sirtl etch time of 5 minutes. 
The dislocation density on the polished side was determined using the standard 
ASTM 9 position measurement technique under a 175X light microscope. Table III 
gives average dislocation densities for the five ASTM points nearest the center 
of the wafer for three typical runs: One wafer polished and inspected on both 

sides (#82) and two wafers (#83 & #84) cycled simultaneously one with polished side 
up and one with the polished side down. 

The quantitative dislocation density fluctuates from run to run but, the wafer 
with the polished side down always has more dislocations than the one with the 
polished side up. The difference in densities from run to run has not yet been 
explained, but the strong correlation between top and bottom dislocation densities 
for a given run is preliminary evidence that gravity influences the creation of 
lattice defects. 


Table III 



AVERAGE DISLOCATION 

DENSITY 

NEAR THE CENTER 

OF THE WAFERS 


Samp le 
Nunib er 
82 

Temperature 

1100" c 

Time 
3 hr 

Polished 

Side 

Both 

Dislocation 

Count 

Top 

Average Dislocation 
Density (per cm^) 
I860 

82 


3 hr 

Both 

Bottom 

4410 

83 

1150“C 

3 hr 

Up 

Top 

435 

84 

1150" C 

3 hr 

Down 

Bottom 

1050 
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III. Theoretical Analysis of the Gravitational Stresses 

Several attempts were made to find a satisfactory means of deteirmining the 
distribution of stress in the wafers supported in the 45® cantilevered position 
under the gravitational bending moments. A computer program was tried which 
uses the finite element method. The results of this technique were found to be 
very inaccurate when a sufficiently small number (200 for a two-dimensional 
analysis) of nodal points were chosen to keep the computing costs at a reasonable 
level. For this reason a closed form calculation was done on a simplified model 
as follows: 

Assume the wafer is a rectangular plate of length, L; thickness, t; and 
width, w: Since the deformation of the wafer is so small the mass (load) distribution 

of the wafer is assumed not to change under gravity and that simultaneous stresses 
may be superimposed. Thus, the stress distribution will be found for a vertically 
standing wafer in a 0.707g field and then for the wafer supported in a horizontal 
cantilevered position in a 0.707g field and the resulting stress distributions 
will be added to obtain the distribution for the wafer supported in the 45® 
cantilevered position. 



VERTICAL WAFER 

let p = mass density of silicon 
X = weight /unit length 
then X 


I0.707g 


.707 p g L w t ^ p g w t 

^ ~ rr 


then cr - x-component of stress for vertical case 

XV 


then 




XV 


wt 
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HORIZONTAL WAFER 




3_x 


+0.707g 


ds El 

where ^ = wafer curvature 

= bending movement 

E = elastic modulus 

I = / y^dA = 1/ 12 w t^ 

also e = 4^ y 
X ds 

where is the strain and y the distance from the 
central plane 


thus cr 


0 = 
X 

Ee = " 

X 

■dd) 

“ Ey 
ds 

= 

\ = 

X(L-x)( 

L-Xv 

2 


: X = 

QgWt 







xh 

pgwt 
2v^ I 

(L - 

x) 

xh " 

6pfi 

(L - 

x) 




45" CANTILEVER 


POSITION 


0=0 + o , 

X XV xh 

^ {(L-x) - ^ (L-x)^} 

^8 =^{(l-x/L) (1-^)2 

t ^ 


} 
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for the surfaces of the wafer y = ± t/2 

'’■.‘bottom) “ ^ ‘ r (1 - 

Notice that the magnitude of the stress at the wafer surface is controlled primarily 
by the ratio L/t. 

3 

for silicon pg - 0,867 Ib/in 

and taking L = 1.5 inch and t = 0.015 inch we get a (5^^ ) = 0.092{fl - —) 

^ x^bottom'^ 

± 300 (1 ^ x/L)^} in pounds/ in^ 

The results are plotted in ELgure II. 

The conclusions to be drawn from this analysis are summarized as follows: 

1. The p^gnit ude of the stresses at the wafer top and bottom surfaces 
are essentially equal, due mostly to the bending moment, and have opposite sign 
(compression for the bottom surface and tension for the top surface) 

2. .The stress is sensitive to the L/t ratio and can be very large for 
tlrin wafers* 

3. The stress of the surface varies slightly less than one order of 
magnitude over the ASTM standard points used in the dislocation count correlation. 
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IV. Correlation of Lattice Defects and Diode Characteristics 

In the course of any study on pro cess- induced defects some attention must 
be paid to their effects on the finished product. This section correlates, 
for a particular diode, the reverse bias breakdown voltage and reverse bias 
leakage current (at 30 volts In darkness) with lattice damage observed by X- 
ray topography. All the processing steps, measurements, and analysis were done 
in the laboratories of the Institute for Solid State Electronics. 

The photograph is a positive enlargement of a Berg-Barrett X-ray topograph 
of the 1.25" wafer used in this study. A recurring pattern is observed with 
a cross at the center of each pattern. The diodes can be identified as the 

rotmd devices located in the lower left comer of each pattern. The bottom of 

the wafer is taken to be the flat near the "J1W8** notation. The upper left 

device is a dual gate FET; the upper ri^t device is an FET; and the lower right 

device is a capacitor. The lattice damage was deliberately produced by very 
rapid thermal cycling of the wafer in a special boat after diffusion and before 
metalization. The defects are most easily seen in the xerox copies near the top 
and sides of the wafer. 

The enclosed wafer facsimiles have blackened squares which indicate the diodes 

that fail to meet a given performance criteria. Most of the "good" diodes have 

a reverse bias breakdown voltage of 80-85 volts with a few as high as 100 volts. 

—9 —10 

Typical leakage currents for the "good" diodes are in the range 10 to 10 
amperes. 

Of the 175 diodes on the wafer, 90 were judged to have some lattice damage 
within the diode spot by careful examination of the original X-ray negative. 



Seventy-nine of the diodes have a reverse bias breakdown voltage of less than 
75 volts. Using 75 volts as the pass-fail criteria the existence of lattice 
defects corresponded to bad devices (or the lack of defects forecast good 
devices) 124 times out of the 175 cases. While this percentage represents a very 
good correlation between defect— device influence, had the defects been produced 
earlier in the process the correlation would probably be even higher. 
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V. Conclusion and Recotnmendations 

First let us review the purpose of this research as it relates to semiconductor 
device fabrication* Thermal processing does induce dislocations through a number 
of well known mechanisms (thermal gradients, oxide Interface stresses, impurity 
gradient stresses, etc.)* However, not all the induced defects can be accounted 
for by these mechanisms and the remainder are generally attributed to "plastic flow"* 

The results on vertically standing wafers conclusively show two significant features of 
the plastic flow phenomenon: 

1* The density of dislocations generated during a given thermal cycle 
is directly related to the duration of the cycle* 

2. The duration of the thermal cycle required to produce a given 
dislocation density is inversly related to the equilibrium 
temperature * 

^v^ile these results by themselves do not prove that gravitationsl stress cause plastic 
flow, it is certainly a plausible explanation. The experiments on wafers cycled in 
the 45^ cantilevered position were done to produce an asymmetrical top to bottom grav- 
itational stress to establish that gravity does play a role in the production of defects 
In examining the results of these 45"^ experiments recall that there is no acceptable 
means of predicting which surface (top or bottom) should have the most dislocations 
under the gravitational bending moments. Thus the fact that there is a contrast in 
the top - bottom dislocation density in wafers thermally cycled at 45*^ is further 
supporting evidence that gravitational stresses contribute to the production of 
defects* The major point is that the effect may not be subtle at alll The theoretical 
analysis of the stress distribution shows that the magnitudes of the stress at the 
top and botto|n surfaces are about the same. The experimental results require the 
conclusion that compressional stress is more likely to produce dislocations* The 
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most difficult feature of the results to understand is the rather small variation 
of the defect density across the face of wafers cycled either vertically or at 
45°. This could be due to the lack of sufficient data for an accurate statistical 
analysis but is more likely due to a more complex distribution of gravitational 
stress than that given by the rather Simple mathematical model. 

Two general conclusions may be drawn from the efforts put into this research: 

1. ) Wafers never reach a true thermodynamic equllbrium condition 

during a thermal cycle. 

2. ) Gravitational stress does play a role in process- induced defect 

generation. 

The most expedient means for firmly and quantitatively establishing the degree 
of influence of gravitational stresses in defect generation is to perform identical 
experiments in other than one-g conditions i.e. centrifuge experiments for greater 
than one-g and orbital experiments for zero-g. It is possible (or even probable) 

’ that defect-free electronic devices could be fabricated in a zero-g environment with 
virtually 100% yields. In view of the small quanities of material involved, zero-g 
processing could have significant economic advantages particularly in the fabrication 


of LSI arrays. 



